Shot noise and reconstruction of the acoustic peak 
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We study the effect of noise in the density field, such as would arise from a finite number density of trac- 
ers, on reconstruction of the acoustic peak within the context of Lagrangian perturbation theory. Recon- 
struction performs better when the density field is determined from denser tracers, but the gains saturate at 
n ~ 10~ 4 (h Mpc -1 ) 3 . For low density tracers it is best to use a large smoothing scale to define the shifts, but 
the optimum is very broad. 
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I. INTRODUCTION 

Baryon acoustic oscillations (BAO) in the baryon-photon 
fluid provide a standard ruler to constrain the expansion of 
the Universe and have become an integral part of current and 
next-generation dark energy experiments These sound 
waves imprint an almost harmonic series of peaks in the power 
spectrum P(k), corresponding to a feature in the correlation 
function £(r)_at ~100 Mpc, with width ~ 10% due to Silk 
damping yj-Q]. Non-linear evolution leads to a damping of 
the oscillations on small scales |&! (and a small shift in their 
positions JM1), 



Pobs(fc) = b 2 e- k ^ 2 / 2 P L {k) + 



(1) 



where we have assumed a scale-independent bias, b, and left 
all broad band and mode-coupling features implicit in the 
The damping of the linear power spectrum (or equivalently the 
smoothing of the correlation function) reduces the contrast of 
the feature and the precision with which the size of ruler may 
be measured and is given by the mean-squared Zel'dovich dis- 
placement of particles, 



(2) 



In 11411 a method was introduced for reducing the damp- 
ing, sharpening the feature in configuration space or restoring 
the higher k oscillations in Fourier space. This procedure was 
studied in lfl5l [Hill using Lagrangian perturbation theory. In 
this brief note we generalize these treatements to show how 
the effects of noise in the density field, arising for example 
from the finite number density of tracers, affects reconstruc- 
tion. We shall concentrate on the broadening of the peak, and 
refer the reader to iT^ fTo"[| for details, discussion and notation. 



field is defined as <5 rcC on = $d — <5s, and to lowest order it 
is equal to the linear density field |[l4T - [l6ll . The non-linear 
damping is however modified from exp[— fc 2 £ 2 /2] to lfl5l,[l6ll 



D{k) = S 2 {k)e 



2S(fc)[l -S(k)]e 



(3) 



with S ss and E^d defined as integrals over the linear power 
spectrum, P L , (see below) and T, 2 d = (1/2) (£ 2 S + Y? dd ). 

If we assume there is a contribution, 6n, from noise we find 
Brecon is unchanged to lowest order. However the damping 
scale is modified. Following 1 1~5| we find 
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d P S 2 {p) [P L (p) + P N (p)} 



(4) 



where Pjy is the power spectrum of Sn and 
1 



3tt 2 



dp [1 - S(p)] z P L (p) + S 2 (p)P N (p) , (5) 



which reduce to the expressions 

of SES 

as Pm — > 0. For 

Poisson shot-noise we expect Pm = &~ 2 n _1 for tracers with 
number density h assuming linear bias b. These equations 
present the generalization of the treament in lfl5l [lql to in- 
clude shot-noise. 



III. RESULTS 

One method to forecast the effect of this noise on cosmo- 
logical parameters constrained by BAO is to replace the Gaus- 
sian damping of Eq. (Q]i with Eq. (O in the computation of the 
Fisher matrix for the acoustic scale s. For example, in spheri- 
cal geometry lfl7ll 



H. RECONSTRUCTION WITH NOISE 

The prescription of [14] begins by smoothing the observed 
density field to filter out high k modes: S(k) — > S(k)8(k). 
We shall take S to be Gaussian of width R. From the 
smoothed field the negative Zel'dovich displacement is com- 
puted s(k) = — i(k/fc 2 )iS(fc)<5(k). Then the objects are 
shifted s to form the "displaced" density field, Sd, and an ini- 
tially spatially uniform grid of particles is also shifted to form 
the "shifted" density field, J s (k). The reconstructed density 



with 
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dP/dlns 
P + n- 1 



P oc D(k) 



sin ks 
ks 



-fc~£g ilk /2 



(6) 



(7) 



where Esuk is the Silk damping scale and • • • refers to terms 
independent of s lfl7ll . The effects of shot-noise show up in 
the increased damping of the higher harmonics of the signal 
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FIG. 1: (Top) The damping scale at z = 0, E e ff, after reconstruction as a function of the number density of tracers, n, assuming 6=1. The 
horizontal dot-dashed line indicates the Silk scale, or the intrinsic width of the acoustic peak, for our cosmology. (Bottom) As above but as a 
function of smoothing scale, R. 



and the increase in the variance per k mode (the denominator 
in Eq.[6]i. 

However, almost as much intuition can be gained by ap- 
proximating D(k) as a Gaussian and asking how the effective 
damping depends on Ppf. To this end we define an "effective" 
S from the value of the damping at fcgd = 0.2 h Mpc - . 

Figure Q] top, shows how £ c ff(-z = 0) depends on n for 
a ACDM model with fi m = 0.25. Note that reconstruc- 
tion improves for higher number density tracers, but the gains 
saturate above approximately 10 -4 ih Mpc -1 ) 3 . For lower 
number densities, it is advantageous to use a larger smooth- 
ing scale to define the shifted field, as expected. For compar- 
isong, without reconstruction the full non-linear smearing at 
z = leads to £ ~ 10 /i -1 Mpc, scaling as the growth factor 
to higher redshift. The horizontal dot-dashed line indicates 
the Silk scale, or the intrinsic width of the acoustic peak, for 
our cosmology — the observed width of the acoustic peak is 
the quadrature sum of Esiik- 

A different view is given in the lower panel of Figure Q] 
which shows how X e ffCz = 0) depends on R for different 



values of n. Note the existence of an "optimal" smoothing 
scale, but that the minimum is extremely broad. 

These results show that, within the context of Lagrangian 
perturbation theory, it is straightforward to understand the ef- 
fects of noise in the density field on the efficacy of reconstruc- 
tion. Reconstruction performs better when the density field 
is determined from denser tracers, but the gains saturate at 
n ~ 10 -4 (ft, Mpc" 1 ) 3 . For low density tracers it is best to 
use a large smoothing scale to define the shifts, but the opti- 
mum is very broad. 
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